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^ \ We discuss lepton flavour violating processes induced in the production and decay 

of heavy right-handed neutrinos at the LHC. Such particles appear in left-right 
' symmetrical extensions of the Standard Model as the messengers of neutrino mass 

generation, and can have masses at the TeV scale. We determine the expected 
D . sensitivity on the right-handed neutrino mixing matrix, as well as on the right- 

handed gauge boson and heavy neutrino masses. By comparing the sensitivity of 
the LHC with that of searches for low energy LFV processes, we identify favourable 
areas of the parameter space to explore the complementarity between LFV at low 
' ' i '■ and high energies. 
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I. INTRODUCTION 



The discovery of neutrino oscillations |lH3j shows that neutrinos are massive |4J and 
that lepton flavour is violated in neutrino propagation. It is natural to expect that the 
violation of this conservation law should show up in other contexts, such as rare lepton 
flavour violating (LFV) decays of muons and taus, e.g. fi~ — > e^7, and possibly also at the 
high energies accessible at the Large Hadron Collider (LHC). In addition to proving that 
flavour is violated in the leptonic sector, oscillation experiments have convincingly shown 
that at least two of the three active neutrinos have a finite mass. However, despite this 
success, oscillation experiments are unable to determine the absolute magnitude of neutrino 
masses. Altogether there are three complementary approaches to probe the absolute scale of 
neutrino mass. Upper limits on the effective electron neutrino mass of ~ 2 eV can be set from 
the analysis of tritium beta-decay experiments 0, 0] . Astronomical observations combined 
with cosmological considerations also allow an upper bound to be set on the sum of the 
three neutrino masses of the order of 0.7 eV, under some assumptions [tJ. The observation 
of neutrinoless double beta decay (Oz^/3/3) would signal the Majorana nature of neutrinos 
and the violation of total lepton number [1,0]. Corresponding searches provide an upper 



bound on the Oz^/3/3 effective Majorana neutrino mass parameter mpp < 300 — 600 meV [10 
Again, it is to be expected that, if Oz//3/3 is observed, the violation of total lepton number 
should also, at some level, take place at high energy accelerators like the LHC. 

Several mechanisms of neutrino mass generation have been suggested in the literature, the 
most prominent example being the seesaw mechanism in which heavy right-handed Majorana 
neutrinos act as "messengers" by generating light Majorana masses for the observed active 
neutrinos through their mixing with the left-handed neutrinos. The Majorana character of 
the active neutrinos can then be connected to a breaking of lepton number symmetry at 
a scale possibly associated with unification [ll- 17] and might also be responsible for the 



baryon asymmetry of the Universe through the leptogenesis mechanism 

Despite its attractiveness, the default type-I seesaw mechanism has important phe- 
nomenological shortcomings: In the standard regime, the right-handed neutrinos have 
masses close to the unification scale and can therefore not be directly produced. In ad- 
dition, the right-handed neutrinos are gauge singlets. This means that even if the masses 
are low enough for them to be produced, the heavy neutrinos only couple with Yukawa 
strength, tightly constrained by the smallness of neutrino masses 1 . This implies that the 



simplest seesaw schemes are difficult to test at the LHC 22 . 

A widely studied alternative of the standard Seesaw scheme with gauge singlet heavy neu- 
trinos is the left-right symmetrical model (LRSM) which extends the electroweak Standard 



Model gauge symmetry SU(2) L ® U(l) to the SU(2) L ® SU(2) iJ ® U(1) B -l group (23j 



26] . Here, right-handed neutrinos are necessary to realize the extended gauge symmetry 
and come as part of an SU(2)# doublet, coupling to the heavy gauge bosons. As a result, 
heavy neutrinos can be produced with gauge coupling strength, with promising discovery 
prospects, given the relatively weak direct experimental bounds on the masses of the extra 
gauge bosons. 

In this paper we address the prospects of probing the mass matrix of the right-handed 
neutrinos. With the existing information on the light neutrino sector and low energy lepton 



1 The mixing of the heavy neutrino singlets is also constrained by precision data such as the rates for lepton 



flavour violating processes 19|-|2l|. 
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(c) 

FIG. 1: Different production processes of heavy right-handed neutrinos with dilepton signatures 
at hadron colliders. 



flavour violation, we re-derive the expected LHC sensitivities on the right-handed gauge 
boson and heavy neutrino masses. Among the possible production processes shown in Fig- 
ure [H we will focus on the resonant production of a right-handed Wr boson displayed in 



Figure 1(a) In addition, we determine the sensitivity of lepton flavour violating processes 
at the LHC to measure the right-handed neutrino masses and their mixing, and compare 
results with low energy LFV searches. The focus in this paper is on lepton flavour violating 
effects rather than lepton number violation, and we will take into account both opposite-sign 
and same-sign lepton events of the process in Figure 1(a) at the LHC. Nevertheless, we will 
briefly comment on dedicated searches of lepton number violation using same-sign event 
signatures and the interplay with Ov/3/3. 

This paper is organized as follows. In Section [III we introduce the minimal LR model, and 
discuss the structure of its charged current couplings. Section HTT1 covers the main observables 
which give current limits on the LR model used. The dilepton signals arising from heavy 
neutrino production are discussed in Section HVl Finally, we present our conclusions in 
Section [V] 



i 



II. LEFT-RIGHT SYMMETRY 

A striking feature of neutrino oscillation data is that the leptonic mixing angles follow a 
pattern substantially different from that which characterizes quarks |27j |. Within the seesaw 
mechanism, as implemented in left-right symmetrical models, neutrino masses are generated 
in a different way from charged fermion masses, as a result of the violation of lepton number. 
Hence, in principle, one may reconcile the large solar and atmospheric angles indicated by 
neutrino oscillation data with the small Kobayashi-Maskawa angles. Many specific models 
based on flavour symmetries have been suggested for this 28[ . In what follows we will simply 



assume a generic lepton flavour structure with the right-handed neutrinos and the charged 
gauge bosons of the left-right model in the range of a few TeV. 



A. Minimal Left-Right Symmetrical Model 

In the minimal SU(2)l <g> SU(2)# (g> XJ(1) b-l model with manifest left-right symmetry, a 
generation of quarks and leptons is assigned to the multiplets (29| 



QL, R =[ u d ) , ^K=VA , 



L,R V / L,R 

with the quantum numbers under SU(2)l (g) SU(2)^ ® U(1) j b-l 

Ql: f^O,^), iJ L : (1,0,- 



2' ' 3/ ' r \2 
Qr: [0,11), 1>R- (0,^-1). (2) 



T 3 J ' r 1 V 2 
The Higgs sector contains a bidoublet 



2' 2 

and two scalar triplets needed to break the right-handed symmetry, 

A ++ 



A " = .f -At". • < 4 > 



with the quantum numbers A L : (1,0,2) and A R : (0,1,2), respectively. This choice is by 
no means unique and there are other multiplet choices that ensure a satisfactory symmetry 
breaking pattern (we will comment on possible alternatives in the conclusion). For our 
analysis we adopt, for concreteness, this more standard choice. The symmetry breaking is 
triggered by the following vacuum expectation values (VEVs) 

= 4) ' {Al ' r) = 6vf o) ' (5) 



where vr of the right triplet breaks SU(2) H ® U(1)b-l to U(l)y and provides masses for 
the new heavy particles. As these new fields have not been observed and due to the strong 
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limits on right-handed currents, Vr should be sufficiently large. On the other hand, the 
VEV Vl of the left triplet contributes to the p parameter, and is therefore experimentally 
constrained to values < 5 GeV (see for example j30|, and references therein). Finally, the 
bidoublet VEVs k\ and k 2 break the Standard Model symmetry and are of the order of the 
electroweak scale. Consequently, the VEVs follow the hierarchy \vj\ <C \ki\ <C \vr\. In the 
following we assume that all VEVs are real and we neglect the possibility of CP violating 
phases. 

The leptonic Yukawa Lagrangian under SXJ(2) L £g> SU(2) fi <g> \](1) B _ L is given by j29[ 



-£ = ^ l YiHr + ^lYiHr + ifl{iY M )&Lil>L + i>l{iY M )A R ^ R + h.c, (6) 

where we suppress gauge invariant field contractions and the summation over the flavour 
indices of the fermion fields and the 3x3 Yukawa matrices Y\, Y 2 arid Ym- 

The Lagrangian (J6]) leads to the 6x6 neutrino mass matrix M v with Dirac mass terms 
arising from Y\ and Y2 as well as Majorana mass terms originating from Ym- The mass 
matrix can be written in block form 

(M L M D \ 

in the basis [y^,v c L ) T . Due to the Pauli principle, M v is complex symmetric. The 3x3 
entries of this matrix are given by 

M L)R = V2Y m v l ,r, M D = ^={Yih + Y 2 k 2 ). (8) 

Assuming that all Yukawa couplings (of a given generation) are of similar magnitude, 
the structure of M v follows the hierarchy of the Higgs VEVs, i.e. M L M D Mr. 
Equivalently, the Yukawa couplings in ([6]) lead to the charged lepton mass matrix Mg = 
■^{Y 2 ki + Yik 2 ). In the following, we will not specify the flavour structure generated by 
the Yukawa couplings. Instead we will phenomenologically parametrize the charged current 
couplings between the heavy neutrinos and the charged leptons. 

The neutrino mass matrix is diagonalized by the 6x6 complex orthogonal mixing matrix 
U v as 

UlM u U u = <\i&g(m Ul ,m V2 ,m V3 ,m Nl1 m, N2 ,m N3 ), (9) 

with the light and heavy neutrino masses m u . and m^., respectively (the transformation can 
be chosen such that the mass eigenvalues are positive and appropriately ordered). 

Before discussing the generic flavour structure of the charged currents generated in the 
LRSM, we will briefly summarize the spectrum of heavy particles relevant to our discussion. 
After both LR and electroweak symmetry breaking and assuming CP invariance, the charged 
gauge boson eigenstates Wl and Wr mix as 

cosCjy s ^w\(Wl\ 
-sinCw cosCwV \WrJ 

Similarly, the neutral gauge bosons Wl, Wl[ and Y combining to Zl and Zr mix to form 
the mass eigenstates 

%\ _ ( cosCz sinCA (Z L S 
Z 2 ) 1-sinCz cosCz/ \Z Rj 




(11) 
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As outlined in Section IHIl the mixing angles (w and (z are highly constrained from elec- 
troweak precision data. Consequently, the mass eigenstates are essentially given by the 
gauge eigenstates, and the heavy components have masses of the order of the right-handed 
symmetry breaking scale, mz R ~ rnw R ~ v r- 

Generally, the neutrino mass eigenstates are mixtures of the weak eigenstates, forming 
three light neutrinos Vi and three heavy neutrinos Ni, cf. Eq. ([9]). Neglecting the flavour 
structure, this leads to the diagonalization in block form 



u\ ( cosCk sin^ 



= ^ bm ^ r L (12) 

^- sin Qu cos Qu J \vr + ^rJ 



with the mixing angle ( u between left- and right-handed neutrinos. For one generation, the 
charged current weak interactions can then be written in terms of the mass eigenstates as 

Jfc = ^= [p + sinC,iV c ) 7 M (1 - 7s)e + ^= sm( w NY(l + 7 5 )e, 

Jfc = ~jj= sinC w ^(l - 7 5 )e + jj= (N - sinC^ c ) Y(l + 7s)e. (13) 

B. Generic Lepton Flavour Structure 

In order to discuss the flavour structure of the charged current interaction we need to 
consider the multi-generation case. The effective 3x6 lepton mixing matrices (Ul^h, 
£ = e, fj,, t, % = 1, . . . , 6 characterizing the charged (and neutral) current weak interactions of 
the mass eigenstate neutrinos in any seesaw model have been fully characterized in Ref. (lEf 
and may be written as 

3 3 

(U L ) H = Kl{U V )ni, {U R ) H = Ke(Uu)n+3,i, (14) 
n=l n=l 

where Q is the 3x3 unitary matrix that diagonalizes the charged lepton mass matrix Mi, 
while U u is the 6x6 unitary matrix that diagonalizes the neutrino mass matrix defined in 
Eq. OH]). For our purposes we take the charged lepton mass matrix in diagonal form 2 so that 
Q — > 1. In this case the overall neutrino matrix U v can be decomposed as 



U » ~ Tri) ~ ttRL ttrr , (15) 



U* L \ _ (U LL U LR 

u R ) - \u RL U RR 

where Ul and Ur relate the left-handed and right-handed neutrino flavour eigenstates vj Jl = 
U£ and VR t = Ng with the mass eigenstates Vi, 

"L,R e = {U L ,R)eM, i = l,...,6, (16) 

and the pieces U LL , U LR , U RL and JJ RR can be calculated numerically or may be obtained 
in seesaw perturbation theory (l6| . 

The charged weak interactions of the light mass eigenstate neutrinos are effectively de- 
scribed by the mixing matrix U LL which is non-unitary, hence the coupling of a given light 



2 This may be automatic in the presence of suitable discrete flavour symmetries as in 
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neutrino to the corresponding charged lepton is decreased with respect to that of the Stan- 
dard Model. This affects the rates for low energy weak decay processes, where the states 
that can be kinematically produced are only the light neutrinos. Similarly, right-handed 



neutrinos would be produced singly in the decays of the Z at the CERN LEP Collider [32 
As a result there are constraints on the strength of the U LR ~ sin mixing matrix elements 
that follow from all these measurements. 

This formalism can be easily adapted to any seesaw model, such as the SU(2)i ® SU(2)# <g> 
U(1)b-l scheme described in Section Hi Al The charged current weak interactions in the left- 
right model can be written in terms of the neutrino mass eigenstates as, 

j£ = i| (u t U^ L + NfU™) T*i L + i?= sin frU RL + N t U^ R ) >ft R , 

= -i| sin (w ipiUt + NiUg") -ft L + ^= (W% R + vtU* L ) -fin, (17) 

in analogy with Eq. (fig)) . Here U LL ,U RR =0(1) while the "small" terms are U LR ,U RL ~ 
MdM^ 1 ~ sin £„, all of which can be obtained within the seesaw perturbative diagonalization 



method developed in Ref. [16J. In the following, we will assume the limit in which all left- 
right mixing terms can be neglected, sin£jy, sin£z, sinCi, <C 1. The applicability of this 
approximation for our calculations will be quantified below. As a result, the only terms 
surviving in (ITT)) are 



J^ R « ^VtiNa^R, (18) 



where we identified U = ?7pmns = U LL and V = JJ RR for notational simplicity. The first term 
describes the mixing of light neutrinos in charged current interactions giving rise to neutrino 
oscillations, whereas the second term is responsible for the LHC process in Figure [U(a) as 
well as all low energy LFV processes as we will discuss below. We will assume manifest 
left-right symmetry of the gauge couplings, i.e. qr = gL, but the general case can always be 
recovered by simple rescaling. 



III. EXPERIMENTAL DATA 



We now turn to the existing experimental constraints on SU(2)l <g> SU(2)# (g) U(1)b-l 
models. Experimental limits on the mass scales and mixing in the minimal LR symmetry 
model come from a variety of sources. For example relevant constraints can be derived from 
the Kl — Kg mass difference, B^B^ oscillations, b quark semileptonic branching ratio and 
decay rate, neutrinoless double beta decay, universality tests, nonleptonic kaon decays, muon 
decay, lepton flavour violating processes and astrophysical constraints from nucleosynthesis 
and SN 1987A 33 35[ . The impact of searches for neutrinoless double beta decay and 



lepton flavour violating processes will be discussed in more detail in Sections IIII Bl and IIII CI 
respectively. Apart from these, the most relevant constraints for our discussions are: A lower 
bound on the Wr mass of mw R > 
from the measurement of the Kl — 
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1.6 TeV 

K$ Kaon mass difference 



37| due to CP violating effects stemming 
but with uncertainties from 



low energy QCD corrections. A more severe limit of mw R > 2.5 TeV was reported in |34l. 135 



s 



and the authors of reference 38j] argue that heavy Higgs masses lighter than 10 TeV are 
disfavoured in the minimal LRSM due to low energy precision data. The minimal LRSM 
is therefore strongly constrained by low energy observations. At the Tevatron, searches for 



Wr — > ev yield a limit of mw R > 1-12 TeV at 95% C.L., assuming SM strength couplings (39 



Only recently, these have been superseded by limits derived from searches at the LHC for Wr 
with decays into ev {mw R > 1-36 TeV) (io| and fiv {mw R > 1-40 TeV) j4l|, again assuming 
SM strength couplings. Using data collected in 2010 and 2011 with an integrated luminosity 
of 240 pb _1 , the CMS collaboration has reported on the search for the production of Wr 
bosons and heavy right-handed neutrinos in the minimal LRSM as in Figure QJa). No excess 
has been observed, thereby excluding a region in the (m\y R , tt^n) parameter space extending 
to (m WR ,m N ) w (1.7 TeV, 0.6 TeV) {42|. Even more recently, the ATLAS collaboration 
has reported on a search for heavy neutrinos and right-handed W bosons via the same 
channel with an integrated luminosity of 2.1 fb~ Again, no excess has been found in 



this search as well, resulting in an excluded parameter region extending to (m\y R ,mN) ~ 
(2.5 TeV, 1.5 TeV). If lighter than « 5 GeV, right-handed neutrinos can be produced on- 
shell in B meson decays, and searches at LHCb provide limits on their coupling strength 
to muons 0, l45f . The bound on the mixing angle between Wr and Wl is of the order 

( w < c(io- 2 ) [HI, M . 

Direct limits on the Zr mass from electroweak precision data, such as lepton universality 



at the Z peak, are of the order 0(1) TeV [47H49| . Within the minimal LRSM one also has 



the theoretical relation rriz R ~ 1.7mw R for qr ~ gi, so that indirect limits via the bounds 
on the Wr mass also yield more stringent constraints on the Zr mass. From the same data, 
the mixing angle between Zr and Zl is constrained to be (z < O(10~ 4 ). 



A. Light Neutrino Oscillation Data 



When combined with reactor and accelerator results, solar and atmospheric neutrino 
experiments [l|-[3| provide firm evidence for neutrino oscillations 0]. Together with reac- 
tor data, the experimental results on solar neutrinos clearly suggest v e — > z/ M . r oscillations 
driven by the mass squared difference Amf 2 = m\ — m\ in the range of the large mixing 
angle solution, while the results on atmospheric neutrinos are interpreted by — > v T os- 
cillations driven by Am^ = m\ — m\ are characterized by a nearly maximal mixing. The 
present global analysis in a three-neutrino framework taking into account the latest results 
on reactor neutrino fluxes gives the following best fit values for the light neutrino squared 
mass differences and the mixing angles #12, #23 and #13 of the neutrino mixing matrix U in 
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the standard parametrization 50j], 



Am 2 21 = (7.62 ± 0.19)10" 5 eV 5 



Am 2 31 



sm 2 9 



12 



sin 2 9 23 



sin 2 6 



13 




-o.ioJ iU 

+0.10N1Q- 

-0.07J iu 



3 eV 2 
3 eV 2 



0.320 

'0.49 
0.53 



+0.015 
•0.017) 

+0.08 
-0.05 
1-0.05 
-0.07 



NH 
IH, 



o.o26i8;88l 
0.02718:881 



NH 
IH. 



NH 
IH. 



(19) 



for normal (NH) and inverse hierarchy (IH), respectively. 



B. Neutrinoless Double Beta Decay 

The most sensitive probe of the absolute mass scale of the neutrinos is neutrinoless double 
(3 decay (0v/3/3). In this process, an atomic nucleus with Z protons decays into a nucleus 
with Z + 2 protons and the same mass number A under the emission of two electrons, 

(A,Z) -> (A,Z + 2) + 2e". (20) 

This process can be engendered through the exchange of a light neutrino connecting two 
V-A weak interactions, as illustrated in Figure Efa). The process in Eq. fl20|) is lepton 
number violating and, in the standard picture of light neutrino exchange, it is only possible 
if the neutrino is identical to its own anti-particle, i.e. if neutrinos are Majorana particles. 
In fact quite generally one may argue that, whatever the underlying mechanism inducing 
Ozv/3/3 decay, its observation implies the Majorana nature of neutrinos jsl, [9[. 

Currently, the best limit on Ou/3/3 decay comes from the search for Ovflfl decay of the 



isotope 76 Ge giving a half-life of Tyi > 1.9 • 10 25 yrs lOj. This results in an upper bound 
on the effective Ov/3/3 Majorana neutrino mass of mpp = | £V U^m Vi \ < 300 — 600 meV, 
depending on the model used to calculate the nuclear matrix element of the process. A 
controversial claim of observation of 0^/3/3 decay in 76 Ge gives a half-life of Ti/ 2 = (O.t 



18.3) -10 25 yrs 5l| and a resulting effective Majorana neutrino mass of mpp = 110 — 560 meV. 
Next generation experiments such as SuperNEMO, GERDA, CUORE, EXO or MAJORANA 
aim to increase the half-life sensitivity by one order of magnitude and will confirm or exclude 
the claimed observation. The planned experiment SuperNEMO allows the measurement of 
Ovfifi decay in several isotopes to the ground and excited states and is able to track the 
trajectories of the emitted electrons and determine their individual energies. In this respect, 
the SuperNEMO experiment has a unique potential to disentangle the possible mechanisms 
for 0vj3f3 decay (52|, 
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In the left-right symmetric model, several such mechanisms can contribute to Oz//3/3 as 
shown in Figure [2] Here, contributions (a)-(d) are due to the exchange of either light or heavy 
neutrinos as well as light and heavy W bosons. All these terms could be described by a single 
Feynman diagram using the mass eigenstates rij (i = 1, 6) and W a (a = 1, 2) of neutrinos 
and W bosons. The separation into the four contributions shown in Figure [2] is traditionally 



10 




"L 



(e) 

FIG. 2: Diagrams contributing to 0v/3/3 decay in left-right symmetry: (a) Light neutrino exchange 
(standard mass mechanism), (b) Heavy neutrino exchange, (c) Neutrino and heavy W exchange 
with Dirac mass helicity flip (A mechanism), (d) Neutrino and light W exchange with Dirac mass 
and W mixing suppression (77 mechanism), (e) Doubly charged Higgs Triplet exchange. 

used as it illustrates the dependence on and suppression with the different LRSM model 
parameters. Diagram 12(a) describes the exchange of massive light neutrinos corresponding 
to the generally considered mass mechanism. Its contribution to CV/3/3 depends on the 
effective neutrino mass mpp = \Yli^ei m ^i\^ an d saturates current experimental bounds if 
the light neutrinos are degenerate with mass scale m vi »s « 0.3 — 0.6 eV. 

Correspondingly, diagram [21(b) describes the exchange of heavy right-handed neutrinos 
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FIG. 3: Diagrams contributing to (J, — > erf. To form complete diagrams, the external photon can 
be attached to any charged particle line. 

and depends on the effective coupling 



If this is the dominant contribution to Ov/3/3, current experimental limits correspond to 
\e N \ < 2x lO" 8 H- 

Diagrams EJjc) and (d) are suppressed by the left-right mixing Mp/M^ ~ Cv ~ \/ m u /mN 
(the latter relation is valid for a dominant type-I seesaw mass mechanism [l^|) between 
light and heavy neutrinos. In our approach we assume that this mixing is small enough 
so that decays of heavy neutrinos via this Yukawa coupling are negligible compared to the 
three-body decays via the SU{2)r gauge coupling. In this case, these contributions to Ovflfl 
are also generally negligible, and we will not discuss them further here. 

Diagram [2|e) describes the contribution from the exchange of a right-handed doubly- 
charged triplet Higgs A^ -3 , with the effective coupling 

A = ^ e *<^<T' (22) 



If dominant, current experimental limits correspond to \e&\ < 8 x 10 54 . 



C. Low Energy Lepton Flavour Violating Processes 

The existence of neutrino oscillations suggests that, at some level, lepton flavour violation 
should also take place in other processes. When taking into account only light neutrinos, 
LFV is strongly suppressed by (AmJ/m^) »s 10~ 50 , due to the GIM mechanism. This 



3 A priori, there is an analogous diagram with a left-handed doubly-charged Higgs, but its contribution is 
always sub-dominant to the standard mass mechanism unless there is a fine-tuning between type-I and 
type-II seesaw contributions to the light neutrino masses. 
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e e 




FIG. 4: Diagrams contributing to fi — > e conversion in nuclei (left) and [i — > eee (right) in left-right 
symmetry. The grey circle represents the effective /i — e— gauge boson vertex with contributions 
from Figure [3 



results in LFV process rates far below any experimental sensitivity which can be safely 
ignored. Within the LRSM, charged lepton flavour violation naturally occurs due to poten- 
tially large flavour violating couplings of the heavy right-handed neutrinos and Higgs scalars 
with charged leptons. Amongst a wide range of possible low energy LFV observables, these 
give rise to observable rates for the processes fi — > ej, \i — > eee and \x — > e conversion in 
nuclei, cf. Figures |3] and HI 

Taking into account contributions from heavy right-handed neutrinos and Higgs scalars, 
the expected branching ratios and conversion rates of the above processes have been cal- 



culated in the LRSM in 55]. In general, these depend on many parameters, but un- 



der the assumption of similar mass scales between the heavy particles in the LRSM, 
m^. ps mw R m m A — ps m A — one can make simple approximations. Such a spectrum 
is naturally expected, as all masses are generated in the breaking of the ri ght- handed sym- 



metry. Under this assumption, the expected branching ratios are given by [55 



R N (fx ->■ e 



T(fj+ ->■ e + 7) 
-» e+vv) 

1.5 x lO" 7 ]^ 
r{n~ + iN -> e~ + %N) 



1 TeV 



(23) 



Br(fi — > eee) 




m, 



a lo£ 



m. 



(24) 



r(/i+ — > e + vv 



1 n I n^w m 
t\h h* I 2 + 

2 \ m . m 



W L 



(25) 



Here, -XVai^Au) ~ (0-8, 1.3, 1.6) is a nucleus-dependent factor whereas g efl and hy describe 
the effective lepton-gauge boson couplings and lepton-Higgs coupling in (quasi-) manifest 
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left-right symmetry, 



^-EU,N 2 > (26) 
hij = y~] V in V jn ( — ^ I , i, j = e, /i, r. (27) 



As shown in [55(, the above approximations are valid if the masses generated in breaking 
the right-handed symmetry are of the same order with 0.2 < mi/rrij < 5 for any pair of 
rrii j = m Nn , mw R , m A£~' 171 a — ■ I n our approach, we keep m WR and m Nn as free parameters 
of the order of 0.5 - 5 TeV, relevant for LHC searches, and consider an order of magnitude 
variation in the heavy Higgs masses of 0.3 < m^--J n ^w R < 3. 

Several properties of the Eqs. fT23|) -( l25|) can then be derived: (i) Both Br(fi — > ej) and 
R N (fi — > e) are proportional to the LFV factor \g ei i\ 2 - In addition, as alog((5 TeV) 2 /m^) = 
0(1), the ratio of their rates gives R N (fi — > e)/Br(fi — > ej) = 0(1), independent of 
the right-handed neutrino mixing matrix V and largely independent of the heavy parti- 
cle spectrum. This consequence of the logarithmic enhancement of the doubly-charged 
Higgs boson contributions to \i — > e conversion is in stark contrast to models where the 
symmetry breaking occurs far above the electroweak scale, such as in supersymmetric see- 
saw models with low tan /3. Here, the Z and photon penguin contributions dominate, and 
R N — > e)/Br(fi — > ej) oc a < C(10 -2 ). (ii) Unless there are cancellations, the LFV cou- 
plings are \g efl \ ~ \h* ee h ep \ and therefore Br(fi — > eee)/R N (fi — > e) = 0(300) (for m A — ps 1 

TeV). 

The above theoretical predictions are to be compared with the current experimental upper 



limits at 90% C.L. 56-58 
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5r cxp (/i -»• erf) < 2.4 ■ 10' 

->• e) < 8.0 • 10" 13 , (28) 
Br cxp (fi eee) < 1.0 ■ 10~ 12 . 

With the current experimental limits roughly of the same order, it follows that the most 
restrictive parameter bounds in the LRSM are derived from the limits on Br(n — > eee). 

As for future developments, the currently running MEG experiment j56[ aims for a sen- 
sitivity of 

Br MEG {fi -> ery) « 10" 13 , (29) 



whereas the COMET and Mu2e experiments both plan to reach a sensitivity of [59|, |60] 

^comet(a< -> e) « 10- 16 . (30) 



IV. DILEPTON SIGNALS AT THE LHC 

In the following we discuss the LHC potential to discover lepton flavour and lepton number 
violating dilepton signals from the production of a heavy right-handed neutrino, 

P + P^W±^£± + N R , (31) 
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FIG. 5: Total cross section of the process pp — > Wr — > Nr£ as a function of the Wr mass, 
for different values of the neutrino Nr mass, calculated with PROTOS 2.1 66(] for the LHC with 
Js = 14 TeV. 



followed by a three-body decay of Nr, as shown in Figure QJa) 



N R ^£f + W R _^£f + 2j. 



(32) 



This process is the main production channel for the right handed neutrinos Nr, where the 
cross section is typically enhanced, compared to other production mechanisms, via on-shell 
Z R production (Figure [TJ^b)) and Wr fusion (Figure [T^c)) [61 64 . 

Since the signal has no missing energy, and since the heavy on-shell Wr subsequently 
decays into an on-shell Nr, backgrounds can be removed and the signal well identified 
by the two Wr and Nr resonances in the respective invariant mass spectra, see Figure [71 
The dilepton signals in left-right symmetric models have already been studied 63M65j . and 

' n. 



bounds on the Wr and Nr masses have been obtained |41|, 

Here we extend the existing analyses and include lepton flavour violation in the right- 
handed neutrino sector to assess the LHC potential to probe the right-handed neutrino 
sector. This allows us not only to determine the Wr and Nr masses, but also to unravel the 
flavour mixing pattern in the heavy Nr sector. These results will be confronted with low 
energy probes of rare LFV processes and neutrinoless double beta decay. 



A. Event Simulation, Topology and Cut Flow 



We generate the partonic signal events qq — > Wr — > £ a N R , followed by the three- 
body decay Nr — > ibjj, with the Triada 1.1 generator 67| of the Monte Carlo package 
PROTOS 2.1 |66|. As a benchmark scenario we choose 



m WR 



2 TeV, 



m NR 



0.5 TeV, 



(33) 



where we assume the equality of the two SU(2) gauge couplings for all calculations, qr = ql, 
and only one of the three neutrinos Nr is lighter than the Wr boson. If other neutrino 
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FIG. 6: Discovery reach of the LHC with y/s = 14 TeV and £ = 30 fb" 1 in the mw R - m N parameter 
plane. The blue solid contours denote constant event rates of the signal process pp —> Wr — > 
££ + 2 jets, after event reconstruction and selection cuts (see Section TlV A 3p . The blue dashed 
contours indicate the discovery reach with S/y/B = 5, and the exclusion region at the 90% C.L. 
The grey solid contour indicates the kinematical threshold mw R = m N for the on-shell decay 
Wr — > £Nr. The grey dashed lines are contours of constant neutrino width T^. Above the red 
solid contours the neutrino decays via an off-shell Wr start to dominate, i.e. Br{N — > £W R ) > 0.5, 
over the SM decay modes N — > £W, uZ, uH [6<3] for a given value of V£ R = sin 2 Cw + sin 2 see 
Eq. (|13p . The red shaded areas are excluded by current LHC searches at CMS [42| (dark shaded) 
and ATLAS 



ij (light shaded). 



states are lighter than Wr, and if their mass difference is sufficiently large, their signals can 
be separated by the resonances in the invariant mass distributions, see Figure [71 Due to 
QCD radiation, the small inherent widths of the neutrinos of some 10 keV get broadened 
up to w 50 GeV. As discussed above, we assume negligible left-right mixing, with (w, (z 
and U LR ~ Q v smaller than 10 -4 [68], such that the SM decays Nr — > W £,Zv,Hv are 
sufficiently suppressed 4 . For the decays of the heavy Wr boson and the neutrino we obtain 
the branching ratios Bv(Nr — > £qq') ~ 83%, for q — u,c, q' — d, s, Bt(Nr — > £tb) w 17%, 
and Bt(Wr — > £Nr) w 9% in our benchmark scenario. 



4 See Ref. 



68( for a discussion of signals with dominating SM decays of Nr. 
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1. Signal Topology 

In Figure we show the mw R dependence of the leading order total production cross 
section cr(pp — > Wr —> ££ + If). For ~ 2.5 TeV and < 1.5 TeV (the exclusion 

reach of current LHC searches), cross sections up to 30 - 50 fb are possible. In Figure [6l 
we show the expected event rates after reconstruction and cuts (as described below) in the 
n^w R - m N R plane. To understand the main kinematic features in the production and decay 
of Wr and Nr, the parameter region can be divided into three different kinematic areas: 

1) The threshold region, with m NR < m WR , is close to where the two-body decay Wr — > 
£Nr is kinematically forbidden, and the signal cross section and thus the LHC sensitivity 
are therefore suppressed. 

2) The jet region, with m^ R <C m WR . Due to the large mass difference, the N R is highly 
boosted. The lepton from the subsequent decay Nr — > t + 2j is not isolated but tends to 
lie inside the cone of the two jets. Although the bare signal cross section is almost constant 
for fixed mw R , the LHC sensitivity in this region is quickly suppressed for < 200 GeV. 

3) The discovery reach zone, with rriN R < rnw R enables ideal decay kinematics with a 
large production cross section and isolated leptons such that the best bounds on the Wr 
mass can be set. 

The blue contours in Figure M give an overview of the LHC reach to probe the Wr and 
mN R masses. As shown by the red contours in Figure El the value of the mixing parameter 
Vlr = sin 2 (w + sin 2 ( u , see Eq. (TlBl) with ql = gR, must be chosen sufficiently small, 
Vlr < 10~ 3 to suppress the SM decays N — > £W, vZ, vH in the relevant parameter space, 



which would otherwise lead to trilepton signatures, see Ref. [68 



2. Background Events 

Since our signal consists of two isolated leptons and at least two high px jets, the major 
SM background sources for opposite-sign dileptons stem from Z + (one or two hard 
jets) and ti production. These backgrounds have total cross sections of the order of 2 x 
10 3 pb [Z + j(j)) and 5 x 10 2 pb (ti), compared to the maximally possible signal cross 
section of ~ 5 x 10~ 2 pb. Other SM backgrounds with subdominant cross sections arise 
from diboson production, pp — > WW, WZ, ZZ, WH, ZH . Note that these SM background 
sources only provide opposite-sign dileptons, whereas same-sign dilepton pairs may originate 
from charge mis-identification, mistakenly reconstructed leptons from jets and SM diboson 
production |69|. which is generally small compared to opposite-sign background. We use the 
code Alpgen [70| to generate the SM Z+jets background, requiring a minimal jet transverse 
momentum of px > 20 GeV to reduce CPU time, and later using proper matching to take 



into account final state radiation. The parton shower Monte Carlo Pythia 6.4 [71[ is used 
to generate the SM backgrounds pp —> ti at leading order. For ti, we take NLO effects into 
account by multiplication with a factor k = 1.6. All hard partonic events are then passed to 
Pythia, to add initial and final state radiation and pile-up, and to perform hadronization. 
Proper matching is applied for initial and final state radiation for the Z+jets Alpgen events. 
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FIG. 7: Distribution of signal events (red) in the invariant masses of the particle combinations 
jj£i@2 ( a ) 5 jj&2 (b), I1I2 (c), as well as the missing transverse momentum (d) in the benchmark 
scenario Eq. ([33]) . at the LHC with 14 TeV and C = 30 fb" 1 . We denote with l^) the leptons with 
the highest (second highest) px- In (c), stacked background events are shown for Z+ jets (light 
green), ti (dark green), and W + Z (blue) production. 



number of jets Nj > 2 

number of isolated leptons Ni = 2 
invariant dilepton mass run > 0.3 TeV 
total invariant mass m Ujj > 1-5 TeV 

TABLE I: Selection cuts used in the LHC event analysis. 
3. Detector Simulation, Event Reconstruction and Selection Cuts 



To simulate a generic LHC detector, we use the fast detector simulation package AcerDET- 
1.0 72] with standard settings, for the simulated signal and background events. As we are 
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TABLE II: Background and signal events after cuts, Table [U for our benchmark scenario, Eq. (|33|) . 
with maximal mixing to electrons and muons, V Ne = Vnh = 1/V2,V Nt = 0, at the LHC with 
14 TeV and C = 30 fb _1 . The dileptons are grouped into pairs with same-sign (SS) and opposite- 
sign (OS) charges. The background event rates for the same-sign signatures were estimated using 
a 5% probability of mis-identifying the charge of one of the leptons. For the signal, we give the 
efficiencies (Eff.) as ratios of events after and before reconstruction and selection cuts in percent. 



interested in lepton flavour violating signals, we generate the event sample using a maxi- 
mal mixing of the active right-handed neutrino to both electrons and muons, but without 
coupling to taus, i.e. we have Vjv e = Vn^ = l/V% and Vn t = 0. For the primary selection, 
we require two AcerDET reconstructed isolated leptons (£ = e, /i) and at least two jets. We 
show the invariant mass distributions of signal and background events before applying cuts 
in Figure [7J 

We then apply our selection cuts as summarized in Table HI to reduce the background [63[ 



64J. Dilepton pairs from Z+jets events are efficiently reduced by requiring a large dilepton 
mass mti > 300 GeV. A cut well above the Z-peak is necessary, due to the long tail in 
the Z — > ii invariant mass distribution, see Figure UJ In addition, a generous cut on the 
total invariant mass mujj > 1.5 TeV reduces tt and Z+jets further, without reducing the 
signal too much, which peaks at mujj ~ mw R -> typically with a width of order 100 GeV 
due to final gluon radiation and smearing. In Figure UJ we also show the signal distribution 
with respect to the missing transverse momentum, but only for illustration, as a cut on 
this variable does not provide a noticeable improvement of the signal over background ratio. 
Nevertheless, appropriate selection criteria with respect to the missing momentum can be 
useful, e.g. in order to decrease possible containment of LNV and LFV signals with SM 
background processes containing light neutrinos. 

For our benchmark scenario, we summarize the background and signal events after cuts 
in Table [Til With these cuts, we obtain signal efficiencies of e ~ 50%, 43%, 34% for the 
dilepton flavour compositions ££' = ee, e/i, fifi, respectively. Note that typically the AcerDET 
reconstruction efficiency is higher for electrons than muons. In addition, the AcerDET algo- 
rithms are not including inefficiencies for the e and \x reconstruction. Additional weighting 
factors of 70 - 90% for each lepton could be applied. On the other hand, we do not include 
a factor k ~ 1.3 for the signal events [73j. In Table [UJ we have also grouped the dileptons 
into pairs with same-sign (SS) and opposite-sign (OS) charges. The background rates of the 
same-sign signatures were estimated using a 5% probability of mis-identifying the charge 



of one of the leptons in the respective opposite-sign signature 63]. No other sources for 
SS background were taken into account. For Wr masses of the order of TeV, the valence 
quarks play a significant role for its production, and thus, reflected by the PDFs, the rate 
for production is larger than that for [63|. Typically the fraction W^/(W^ + Wr) 
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of produced W£ in pp collisions at 14 TeV changes from »s 70% to ~ 95% for m WR in- 
creasing from 1 TeV to 10 TeV j63|, We thus obtain for the reconstructed signal events 
N(e+e-) : N(e + e + ) : N(e~e-) « 7V(/x+/i~) : N(fj, + (j+) : iV(/x"//-) « V(e + /i~ + e"/i + ) : 
iV(e + /i + ) : N(e~ p~) ~ 4 : 3 : 1, in our benchmark scenario. 



^. Selection Cut Optimization 

For the purpose of our study, we do not optimize the cuts for each LRSM parameter point 
defined by mw R an d ttin r - Instead, we always use the cuts as given in Tabled which were 
chosen to highlight the discovery reach for the LHC The cuts could be adapted by requiring 
higher transverse momenta of the leptons and jets, pt > 100 GeV as well as selecting on 
missing p™ lss < 50 GeV [68[. Further, one can apply mass- window cuts to isolate the Wr 
and Nr mass peaks in the m^j and Vfiua invariant mass distributions, respectively [63|, 64]. 
In particular, different m^j mass-window cuts might be necessary to disentangle different 
signal contributions, if more than one Nr is kinematically accessible. Also, our cut on the 
total invariant mass mujj > 1.5 TeV is chosen for heavy bosons m-w R > 2 TeV, and for 
smaller masses one loses signal events, which can be observed in the indentation on the left 
side of the 10 3 events contour in Figure [6j 

The cuts might also be improved in the kinematically suppressed parameter regions, i.e. 
for mN R < mw R (threshold region) and mN R <C mw R (jet region). To enhance the sensitivity 
in the threshold region, same-sign lepton pairs can be selected, or the cuts can be raised 
for the px of the jets and for the dilepton invariant mass m« [63|. In the jet region, where 
the lepton from the Nr decay tends to lie inside the jet cones, one could search for events 
with one high p? isolated lepton and one high p? hadronic jet with a large electromagnetic 
component, and matching a high p? track in the inner detector. This typically increases the 



sensitivity in those regions by a few percent |63l . [6 



For same-sign dileptons, the kinematic cuts can be relaxed, since SM background events 
will dominantly produce opposite-sign dileptons. As mentioned above, same-sign lepton 
background events would originate from charge mis-identification, mistakenly reconstructed 
leptons from jets and SM diboson production j69 ! ']. The dilepton pairs could also be grouped 
according to their flavour configurations; for example, lepton pairs with different flavours 
only originate from ti production, reducible by vetoing 6-jets. Since the b-tagging efficiencies 
at the LHC are of the order of 60% to 70%, the ti background should be reduced by an order 
of magnitude. The signal loss will be small, since we typically have Bt(Nr — > Ibt) ~ 20%, 
and Br(t — > btv) = 20%, £ = e, /x. However, as pointed out earlier, to estimate the discovery 
potential of the LHC for lepton flavour violating signals, we will base the following analysis 
on the basic conservative kinematic cuts as given in Table [B 



B. Lepton Flavour Violation 

In this section we discuss the prospects of observing lepton flavour violating production 
and decay of right-handed Wr bosons and heavy neutrinos in the left-right symmetric model 
described in Section III A\ due to a mixing among the right-handed neutrinos. We will focus 
on a mixing of right-handed neutrinos with electrons and muons, i.e. we are mostly interested 
in the mixing matrix elements V^ ie and V^ ifJL . In turn, we will assume the case of flavour 
mixing with the unitarity constraints V^. e + = 1 (i = 1,2), Vjv ie Vjv 2 e + VN lfJ ,VN 2 ^ — 0) 
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FIG. 8: The solid blue contours denote the event rates for the processes pp — > Wr — > e^/i^ + 2 jets 
(left) and pp — > Wr — > e ± / u ± ' =F + 2 jets (right), as a function of right-handed W boson mass mw R , 
and the right-handed neutrino mass mjy, at the LHC with 14 TeV and C = 30 fb" 1 . The grey 
contour corresponds to the kinematical threshold mw R = m^. The blue dashed contours give the 
discovery and exclusion reach with S/^fB = 5 and 1.64 (90% C.L.), respectively. The processes 
are calculated for maximal unitary coupling of the right-handed neutrino to e and fx only with 
V^e = = 1/V2. The red shaded areas are excluded by current LHC searches at CMS 42 1 
(dark shaded) and ATLAS 0] (light shaded). 



with either one or two heavy neutrinos light enough to be produced at the LHC 5 . In the first 
case, we will also discuss the more general scenario of a possible mixing to taus by extending 
the unitarity relation as V^ ie + V$ + V^ lT = 1. Here, we will not take into account taus in 
the final state 6 , but a non-zero value of Vjv ir will reduce the mixing to electrons and muons 
accordingly. 



1. Single Neutrino Exchange 

To simplify our discussion we first consider one right-handed neutrino in the intermediate 
state. That is, either only one right-handed neutrino is light enough to be produced, m Nl < 
mw R < m N 2 ai 01 if more than one neutrino is below the threshold, the mass difference 



5 We always assume real mixing matrix elements, i.e. we neglect any possible CP violating phases in the 
right-handed neutrino sector. 

6 This is a highly interesting possibility in its own right as the reconstruction efficiencies in the leptonic 
channels r — > lv(V r with I = e, /i are only reduced by about a third of those for electrons and muons, 
simply due to the leptonic branching ratios of the r 7J| . This could make it possible to detect r — \x and 
t — e flavour violation at the LHC. 
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between the right-handed neutrinos is sufficiently large, Am^^ > 100 GeV, such that the 
neutrino resonances can be individually reconstructed. In Section IIVB2[ we then discuss 
right-handed neutrinos with smaller mass differences. 

The couplings of the right-handed neutrino Ni = N to electrons and muons are given by 
the elements Vjv e and Vjv> of the right-handed neutrino mixing matrix. After reconstruction 
and cuts, the number of opposite-sign (OS) dilepton events from the processes pp — > Wr — > 
e + e~(/i + /i~, e ± / u =F ) + 2 jets, are given by 

N(e + e-) = e°J xCx a e ° e s (m WR ,m N ) x |V^ e | 4 , 

NQj+h-) = e™x£x a° e s (m WR , m N ) x \V N ,\\ (34) 
iV(e ± /i T ) = e°* x Cx a° s (m WR ,m N ) x | VAr e | 2 1 Vat^ | 2 , 

with the integrated luminosity £, the cross section cr° e s (mw R , m N) for e + e~ production 
with Vjve — 1 ( n ° flavour mixing), and the experimental efficiencies e° e s , e°^, due to 
reconstruction and cuts, for the different dilepton flavour combinations. Similar relations 
are obtained for same-sign dilepton pairs, with the efficiencies e^, ejf , eft and the SS cross 
section (m\v R ,mpi). The numerical values of the efficiencies are summarized in Table HTl 
obtained in our simulation of the benchmark scenario mw R = 2 TeV, = 0.5 TeV with 
the cuts shown in Table [B 

To determine the reach of the LHC to probe lepton flavour violation, we show in Figure [8] 
the event rate for the signature pp — > Wr — > efi + 2 jets at the LHC with 14 TeV and an 
integrated luminosity of C = 30 fb _1 for maximal mixing V^ e = V^^ = l/y/2. In Figured] 
(left) we show the sum of events of opposite-sign leptons, iV(e + /!~)+iV(e~/! + ), whereas in the 
right plot the charges are summed over as well, N(e + fi~) + N(e~ n + ) + iV(e + /i + ) + N(e~fi~). 
The blue dashed contours show the discovery reach {S/\fB = 5, practically coinciding 
with the 100 events contour) and the exclusion reach (90% C.L. excess over background). 
As discussed previously, we do not optimize the cuts for each individual parameter point, 
which could extend the reach somewhat. Because of the simple relations in Eq. (I3~4"|) . the 
LHC reach is comparable for lepton flavour conserving and maximally violating signals, c.f. 
Figures [6] and Effectively, the rate is halved for maximal mixing, and the shape of the 
regions follows from the discussion in Section IIV1 It should therefore be possible to probe 
LFV for Wr masses up to 3-3.5 TeV and right-handed neutrino masses up to 1.5-2 TeV. For 



exclusion region with C = 2.1 fb _1 US 



I up 

comparison, we show the CMS exclusion region with £ = 0.24 fb _i |42|, and the ATLAS 



To verify the LHC sensitivity to the flavour couplings V^ e and V^^, we fix the masses 
mw R = 2.5 TeV, = 0.5 TeV, and show in Figure M the event rates as a function of 
the electron coupling V^ e for OS and OS+SS signatures, where we assume unitarity among 
the couplings Vjv e and Vn^, i.e. V^ e + = 1. The dependence of the event rates on the 
couplings is again very simple and follows Eqs. fl34l) . If the lepton flavour violating signature 
e-/i is observable, at least one of the flavour conserving signatures e-e or can be observed 
as well. If V^, > 1/2, the neutrino couples more strongly to electrons than to muons, and 
from Eqs. fl34"l) we can determine 



V 2 
> t 



Nc 



l6f/ iy(e+/i- + e-/i + ) 
+ 2eOS N(e+e~) 



(35) 



for opposite-sign signatures. The coupling parameter | Vjv e | can therefore be determined 
independently of the absolute normalization of the production cross section (e.g. due to 
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FIG. 9: Event rates for the processes pp —> Wr — > ee^fifj,, efi) + 2 jets at the LHC with 14 TeV 
and £ = 30 fb _1 as a function of the right-handed neutrino-electron coupling V^ e = 1 — Vjy for 
a right-handed W boson mass mw R = 2.5 TeV and right-handed neutrino mass mjy = 0.5 TeV. 
The rates are calculated using only the opposite charge sign event sample (left) and both opposite 
and same charge sign sample (right) of the two leptons. The horizontal dashed lines show the 
background to each correspondingly colored signal. 



possible splitting of the gauge couplings, Qr ^ Ql effects) and systematic uncertainties are 
reduced. If V^ e < 1/2, we find an analogous expression, replacing N(e + e~) with N(fi + fi~) 
in the denominator of Eq. (1351) . For V^ e ~ 1/2, Eq. (155]) can be used both with N(e + e~) 
and N(fi + fi~), yielding two independent measurements for V$ e . 

This can be generalized by dropping the unitarity condition and keeping the couplings 
V/Ve and Vn^ independent from each other. This approach is used in Figure [TU1 where we 
show the excluded region in the V^ e and V/v M parameter plane at 90% C.L., using a x 2 
analysis of the ee, e/i and \x\x events with both OS and OS+SS signatures, assuming no 
excess above the background is observed. In the given scenario with mw R = 2.5 TeV and 
771 7v = 0.5 TeV, couplings of the order of V^ e ^ ~ 0.5 can be excluded. For illustration, 
we also show uncertainty contours for hypothetical signals corresponding to four different 
choices for (Vjv e , Vjy^): (a) (l/y 7 ^, l/y/2), corresponding to maximal unitary mixing; (b) 
(0.347,0.938), representing the unitary scenario with the minimal value of Vjve that can 
be distinguished from zero; (c) (0.937,0.348), representing the unitary scenario with the 
minimal value of that can be distinguished from zero; (d) (0.5, 0.5), representing a non- 
unitary scenario close to the exclusion limit. In Figure UJJ\ the errors on the event rates are 
assumed to be dominated by their statistical uncertainties. 

In order to derive the sensitivity as a function of mw R and m/v, we show in Figure [TT] 
the minimal coupling V^ e which can be observed at 5a, for both OS and OS+SS leptons. 
The outermost contour in the plots corresponds to a signal at 90% C.L. for maximal mixing 
V£ e = = 1/2. This exclusion contour and the 5a discovery contour therefore correspond 
to the dashed contours in Figure 13 respectively. The parameter region with the largest cross 
section around (mw R , win) — (1-6, 0.3) TeV is already excluded by LHC searches, and flavour 
violating right-handed neutrino-lepton couplings down to V^ e ^ ~ 10 _1 can potentially be 

probed at the LHC with 14 TeV and C = 30 fb _1 . 
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FIG. 10: Excluded regions in the Vj\r e and Vjv> parameter plane at 90% C.L. using OS (light red 
shaded with solid boundary) and OS+SS (dark red with dashed boundary) event samples, assuming 
no excess above background. Also shown are the la and 5a uncertainty solid blue contours for 
hypothetical signals corresponding to (V Ne ,V Nll ) = (l/v/2, 1/V2), (0.347,0.938), (0.937,0.348), 
(0.5, 0.5) using OS+SS samples. The statistical analysis is based on the events pp — > Wr — > 
ee(fi/j,, e/x) + 2 jets at the LHC with 14 TeV and C = 30 fb" 1 , with a Wr boson mass of mw R = 
2.5 TeV and right-handed neutrino mass of rriN = 0.5 TeV. The errors on the event rates are 
assumed to be dominated by their statistical uncertainty. 



2. Two Neutrino Exchange 



In general, all heavy neutrinos with < mw R that couple to electrons and/or muons 
will contribute to our dilepton signatures and have to be taken into account. When summing 
over the heavy right-handed neutrinos in the intermediate state, this may lead to interference 
effects between the different contributions. In fact in the limit of degenerate heavy neutrinos, 
Am-j = m 2 Ni —m 2 N . — > 0, all lepton flavour violating signals will suffer a GIM-like suppression, 
analogous to low energy LFV processes as described in Section IIII CI As a crucial difference 
to the radiative rare decays, the neutrinos are produced on-shell at the LHC with a short 
decay length. This leads to a decoherence of the right-handed neutrino oscillation, and the 
suppression is proportional to Am 2 j / (ttinT n) , rather than Am 2 Jm 2 N . This follows from the 
well justified narrow width approximation for the product of the neutrino propagators JV$ 
and Nj in the squared matrix element 
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FIG. 11: Sensitivity to the coupling as function of mw R and m^ R at the LHC with 14 TeV and 
C = 30 fb _1 using OS (left) and OS+SS leptons (right). The solid contours indicate a discovery 
sensitivity at 5a and the outermost contour corresponds to an exclusion at 90% C.L. for maximal 
mixing Vj^ e = 1/2. The grey contour corresponds to the kinematical threshold mw R = n^N- 



The red shaded areas are excluded by current LHC searches at CMS 42| (dark shaded) and 
ATLAS (light shaded). 



with 

/ Am 2 - \ 1 1 
Cij = \ 1 + i 2(mT) ) ' Am ^' = m ^ ~ = 2^ mNiTNi + mN ^ N ^ ^ 



For large mass splittings, Arn-j ^> (mr).., the factors approach Sy, i.e. the neutrino 
interference is suppressed and their contributions add up incoherently in the squared matrix 
element. For small mass splittings, Am| <C (ml 1 ).., due to Eq. fl36|) and the unitarity of 
the mixing matrix V in the vertex V^eNi-i-Wpt, the flavour violating process pp — > Wr — > 
£\£ 2 + 2 jets, £i ^ £ 2 , is suppressed as Am? / (mT)^. When two neutrinos iV li2 couple to e 
and /x, with a unitarity matrix described by the mixing angle <fr, 

v= (cos^ sin<A 
\— sin0 cos0/ 

the total event rate for the process — >■ — >■ e ± /i =F + 2 jets is then 

1 (Am 2 ) 2 

V(eV) ~ e ^C[a(m WR ,m Nl ) + a(m WR ,m N2 )] x -sin 2 (20) 2 r , 2 • ( 39 ) 

Here, a (m^ fl , mjvj is the cross section of the lepton flavour conserving process pp — > e + e~ + 
2 jets with exchange of one neutrino of mass m^ coupling only to electrons. This equation 
describes the event rate for all mass splitting regimes to a very good approximation: 
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1. If the neutrino mass difference is larger than the experimental width of the neutrino 
resonances in the invariant mass distribution rri£ 2 jj, Arriij > T e ^ , it is in principle 
possible to reconstruct each of the resonances m Nl2 , and to determine the couplings 
^/Vie' ^TViAf ^jv 2 e' ^N 2 fi °f both neutrinos, as described in Section HVB 1[ The individual 
event rate for each resonance is then 

^(e ± /i T ) ~ e™£(T(mw R ,rn Ni ) x ^sin 2 (20), 

which is independent of the neutrino mass splitting. 

2. If r^y. <C Arriij *C T^, it is not possible to resolve the individual neutrino contribu- 
tions which overlap and form a single resonance at « « m^ 2 . Nevertheless, 
the LFV rate is still unsuppressed by the mass splitting and the total event rate in 
the single resonance is given by 

N(e ± ^) « e^Ca(m WR ,mjf) x ^sin 2 (20). 

3. If Arriij "C TjVj, the neutrinos form a single resonance at and the total event rate 
is suppressed with Amy = m N . — m Nj as 

/Am \ 2 

iV(eV) » 6^ s £a(m^,m^) x sin 2 (20) (— , (r> = (r^ +r^)/2. 

Due to the small neutrino three-body decay width, IV oc g\rri° N j rri^ < 10~ 2 GeV (cf. 
Figure [6]), the regimes 1) and 2) apply for a wide range of neutrino mass splittings. In 
these regimes, the oscillations between right-handed neutrinos decohere. Consequently, the 
LFV process rates at the LHC are not suppressed and are independent of the mass splitting 
(unless the individual resonances can be resolved kinematically) , and can be probed at the 
LHC for mass differences as small as Am^ « « 10~ 2 — 10 -6 GeV. This is in stark 
contrast to low energy rare LFV processes which, as described in Section IIII C\ experience 
a right-handed GIM-like suppression as Amfy/m^ . For example, for two flavour e — \x 
mixing, the branching ratio of fi — > e'y, Eq. (1231) . can be written in our parametrization as 

q s /Am 2 \ 2 /2 TeV\ 4 , s 

fir( M ->e 7 )^2xlO- 9 sin 2 (20) -— , (40) 

V m w R J \ m w R J 

with similar results for the other rare processes. For the unsuppressed case with large 
mixing m tt/A and large mass difference Am 2 2 > rn^^ the current limits on this decay 
already put severe constraints on the scale rriw R - As discussed in Section UlI C[ the // — > e 
conversion in nuclei and fi — > eee are expected to restrict the parameter space even more. 
On the other hand, even for modest neutrino mass splittings, e.g. m^ 2 — « 50 GeV in 
our benchmark scenario with — 500 GeV and rriw R = 2000 GeV, the branching ratio 
is already suppressed by a factor « 10 -4 to Br(fi — > ej) ~ 4 ■ 10~ 13 , below the current 
experimental limit. The LFV event rate at the LHC would not be affected by a mass 
splitting of this size, and it could even be possible to reconstruct the individual neutrino 
mass resonances. Measurements of LFV processes at the LHC and low energy rare LFV 
decays therefore provide highly complementary information on the mass spectrum and the 
flavour couplings of the LRSM. 
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m WR [TeV] 

FIG. 12: Comparison of LFV event rates at the LHC and in low energy rare decays. The solid 
blue contours give the number of OS+SS events for the LFV signature e =t /x =l: ' =F + 2j at the LHC 
with 14 TeV and C = 30 fb _1 . The dashed contours define the parameter region with signals at 
5(7 and 90%. The shaded red areas denote the parameter regions excluded by current low energy 
LFV limits, whereas the red contours show the expected sensitivity of planned experiments. The 
processes have been calculated using maximal unitary mixing between two heavy neutrinos coupling 
only to e and fj, with a 1% mass splitting, i.e. <p = ir/4 (cf. Eq. [38]) . (mjv 2 ~~ m 7Vi)/ m A r i = 0.01. 
The spectrum of the doubly charged Higgs bosons is given by m A — = mw R - The grey contour 
corresponds to the kinematical threshold m\v R = m^. 

To explore this complementarity, we first compare the sensitivity of \i— e LFV processes on 
the masses mw R and tun for maximal mixing <fi = n/A and a fixed 1% neutrino mass splitting, 
(mjv 2 — ttiatJ / 'ttitv = 10~ 2 . This is shown in Figures fT2l and [T3l for different heavy Higgs boson 
mass spectra: i m^-- R l mw R — 0.3,1,3, mapping out the naturally expected spectral range 
of scales in the LRSM. As discussed in Section fill CI the lower the Higgs boson masses are, 
the more constraining the processes \i — e conversion and \x —> eee become. The current 
limits on the rare processes already strongly constrain the parameter space, with \x —> eee 
providing the most stringent bound. Even for m A — /t^Wr — 3, the parameter space with 

rriN > rn\v R i which is inaccessible at the LHC, is almost ruled out. For m A — /mw R — 0.3 
the parameter space that could be probed at the LHC is almost ruled out, and only the 
region with rather low neutrino masses, m^r ~ 200 — 600 GeV, is still allowed. In this 
mass regime, the expected COMET sensitivity on R Al (fi — >■ e) = 10~ 16 will fully probe the 



27 




m WR [TeV] m WR [TeV] 

FIG. 13: As Figure [121 but with a doubly charged Higgs boson mass spectrum m A — = 0.3 x mw R 
(left) and m A — = 3 x mw R (right). 

LHC accessible parameter space. Generically, there is a high potential that low energy LFV 
processes can probe the parameter space that can be also tested at the LHC, allowing for a 
highly detailed view of the LRSM mass spectrum and flavour mixing properties. 

To explore the complementarity in probing different neutrino mass splittings and flavour 
mixing, we compare the sensitivity of low energy and LHC processes on sin 2 (20) and 
AmN/rriN in Figure HH Here, the LRSM mass spectrum is fixed to = 0.5 TeV, 
mw R = 2.5 TeV and m A — = 0.66 TeV, corresponding to a light Higgs sector. As 

discussed above, the LFV process rate at the LHC is independent of the neutrino mass 
splitting until it becomes comparable and smaller than the heavy neutrino decay width at 
r^v/mTv ~ 5- 10~ 8 . It is therefore possible to probe such tiny mass differences at the LHC for 
mixing angles sin 2 (20) > 10 _1 in this scenario. On the other hand, the low energy processes 
exhibit the typical dependence oc sin 2 (20)(Am 2 v ) 2 , and may only probe mass splittings as 
low as Am7v/m7v ~ 10~ 3 . 

C. Lepton Number Violation 

So far we have only considered OS and OS+SS event signatures, but not the SS sample 
independently, as the significance of such a signal depends crucially on a proper treatment 
of the same-sign lepton background arising from charge mis-identification, mistakenly recon- 
structed hard leptons from jets and diboson production. Such a treatment goes beyond the 
scope of this work. Same-sign lepton events are of course a crucial consequence of the Ma- 
jorana nature of heavy neutrinos in left-right symmetric models, and the associated lepton 
number violation. In the LRSM considered here, the heavy neutrinos acquire their Majorana 
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FIG. 14: As Figure [T21 but showing the dependence on the heavy neutrino mixing angle parameter 
sin 2 2(f> and the mass splitting (mjv 2 — n^N^/^Ni- The mass scales are given by mw R = 2.5 TeV 
and mjvi = 0.5 TeV, with a light heavy Higgs sector, m A — = 0.3 x mw R - The dashed horizontal 
line denotes the value for the heavy neutrino width Tn/ttin. 



masses through the breaking of both the left-right and lepton number symmetry at the high 
scale Vr, much larger than the neutrino Dirac mass terms, see Eq. (JHJ). 

As a result, processes mediated by the heavy neutrinos exhibit maximal lepton number 
violation, with equal probabilities for the heavy neutrinos decaying into positively and nega- 
tively charged leptons. Such a scenario is technically natural to describe both light neutrino 
masses « 0.1 eV and heavy neutrino masses ~ 1 TeV, if one assumes ad hoc small Dirac 
couplings between light and heavy neutrinos. On the other hand, accounting for the light 
neutrino masses within the seesaw mechanism with heavy neutrinos close to the electroweak 
scale would require the breaking of lepton number at a much lower scale /^lnv 5 as in the 



inverse or linear seesaw schemes [76 H80I. o r by invoking a proper flavour symmetry among 



the lepton Yukawa couplings, see e.g. |8ll4S4j. These scenarios however lead naturally to 



quasi-Dirac heavy neutrinos where all lepton number violating processes are suppressed by 
the small mass splitting of heavy neutrinos, /iLNv/ m 7v- From the model building viewpoint 
it is therefore not straightforward to predict the rate of same-sign lepton events for a given 
heavy neutrino mass, whereas the rate of events with opposite-sign but different flavours are 
generally independent of the mechanism of lepton number violation. This question will be 
addressed in an upcoming analysis. 

Despite the above caveats, for illustration we compare the sensitivity of observing lepton 
number violation at the LHC with that of Of/3/3 experiments in the LRSM. As pointed out in 
when staying within the LRSM, it is possible to explore the interplay between LFV 
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FIG. 15: Comparison of LNV event rates at the LHC and in 0^/3/3 experiments. The solid blue con- 
tours give the number of same-sign events for the LNV signature e^e^ + 2j at the LHC with 14 TeV 
and C = 30 fb" 1 . The dashed blue contours define the parameter regions with a significance at 
5<7 and 90%, respectively, where the background is estimated using a 5% charge mis-identification 
probability (see discussion of Table UTJ) . The shaded green area denotes the parameter space ex- 
cluded by 0v(3(3 at T ^ w 2 x 10 25 yrs, assuming dominant doubly-charged Higgs or heavy 
neutrino exchange. Correspondingly, the green dashed contours show the sensitivity of future 
Ouf3f3 experiments at T ^^ f» 2 x 10 26 yrs. The red shaded area is excluded by current LHC 
searches at ATLAS 



and LNV further. The LHC analysis in this section is the same as described in Section HV A[ 
using identical cuts and reconstruction criteria. As illustrated in Table [TTJ we estimate the 
same-sign dilepton background by assuming a 5% probability of mis-identifying the charge 
of a lepton of an OS event. 

In Figure PL31 we show the event rate of the same-sign signature e ± e ± + 2j at the LHC 
with 14 TeV and £ = 30 fb" 1 . For comparison, the shaded green regions and green dashed 
contours represent the current limit from T "^ 13 > 2 x 10 25 yrs (Heidelberg- Moscow) and the 
future sensitivity T 0l//313 m 2 x 10 26 yrs (improvement by one order of magnitude) of Qv/3/3 
experiments. Here, we assume that neutrinoless double beta decay is either dominated by 
heavy neutrino exchange (Figure [2(b) and described by the effective coupling in Eq. (I2ip ) or 
by Higgs triplet exchange (Figure [2(e) and described by the effective coupling in Eq. (I2"2"j) ). 
In the latter light doubly charged Higgs boson mass spectrum is assumed, with 

m A — = 0.33 x mw R - in both cases any other contribution is assumed to be zero. As 
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discussed in Section IIII B[ this is naturally the case for the left-right mixing contributions 
Figure [2](c) and (d), though the contribution from light neutrino exchange (Figure 12(a)) will 
always be present. Hence the sensitivities shown here would correspond to an hierarchical 
light neutrino mass spectrum, with a small effective mass mpp. Due to the uncertainties 
described above and because of the additional dependence on the doubly charged Higgs mass 
and the heavy neutrino mixing matrix (here we neglect lepton flavour violation, i.e. we use 
Vjve = 1, = Vn t = 0), the sensitivities shown here are for illustration only. Nevertheless, 
it is interesting that the area currently probed by the LHC will also be tested in upcoming 
neutrinoless double beta decay experiments. 



V. CONCLUSIONS 

Light neutrino masses naturally arise in left-right symmetric seesaw extensions of the 
Standard Model, as required in order to account for current neutrino oscillations data. 
Lepton flavour violating effects involving charged leptons are also naturally expected in such 
scenarios if the masses of the heavy right-handed neutrinos present in left-right symmetric 
models are of the order of 1 - 10 TeV. Here we have considered lepton flavour violating 
processes induced in the production and decay of heavy right-handed neutrinos and the 
resulting signatures at the LHC. Either through the assumption of small Yukawa couplings or 
suitable model constructions implementing a low-scale seesaw, for example inverse or linear 
seesaw, right-handed neutrinos can have masses of order TeV, and are hence accessible at 
the LHC. For this case we have derived the expected LHC sensitivities on the right-handed 
gauge boson and heavy neutrino masses, as well as the LFV couplings of the heavy neutrinos 
to charged leptons and compared the collider results with existing bounds from low energy 
LFV rare decays. Our discussion was mainly devoted to the first two leptonic flavours, due 
to their cleaner detection prospects, though extension to the tau flavour will be important, 



in view of good tau detection efficiencies in decays of heavy particles at the LHC [74 . 

Because right-handed neutrinos can be produced at gauge coupling strength in left-right 
symmetric models, the LHC has the potential to discover right-handed Wr bosons up to 
m w R ~ 3 — 4 TeV and heavy neutrinos up to m^ R « 2 — 2.5 TeV, at 14 TeV and C = 30 fb" 1 . 
No signal has been found so far, and current bounds from LHC searches are already stringent, 
with the excluded area extending to (m WR ,m NR ) « (2.5,1.5) TeV. Outside these limits, 
LFV couplings, described by the heavy neutrino mixing matrix V entering the right-handed 
charged current interaction ^VmiN r ^I r , as low as V^i ~ 0.3 can still be probed at the 
LHC. We have also explored the complementarity of such searches with LFV probes at low 
energies, namely n — > ej, ji — > 3e and /i — > e conversion in nuclei. If the mass splitting of 
the heavy right-handed neutrinos is large, Am 2 N /m 2 VR > 1, these processes already heavily 
constrain the presence of LFV in the right-handed neutrino sector. On the other hand, if 
the mass splittings are small, Am 2 v /m 2 VR < 0.01, the low energy LFV processes are GIM 
suppressed whereas LFV can still be observed at the LHC through the resonant production of 
right-handed neutrinos. This analysis provides an example of the general complementarity 
between LFV searches at the LHC and at low energies. Being based on high intensity 
experiments, low energy probes such as /i — > cy have a further reach to higher scales of new 
physics as well as smaller LFV couplings and mixing angles. On the other hand, collider 
searches are limited by the available energy and luminosity, but have the potential to probe 
individual particles and couplings. 

In this work, our focus has been on lepton flavour violating effects and the potential 
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FIG. 16: Comparison of LFV and LNV event rates at the LHC and in low energy probes. The solid 
blue contours define the parameter region with signals of the LFV OS+SS process pp —> Wr — > 
N R l -»• e/i + 2 jets at the LHC at ha and 90% (14 TeV, C = 30 fb _1 ). The solid gray contour 
corresponds to the kinematical threshold mw R = m N- Overlayed are the current and expected 
future sensitivities of low energy LFV processes and Oz^/3/3 (mediated by heavy neutrinos), as 
denoted in the plot. All processes were calculated assuming maximal flavour mixing of two heavy 
neutrinos N± and N2 to electrons and muons with a mass difference rn/v 2 — mN-i = O.Olmjv- The 
red shaded area in the lower left corner is excluded by current LHC searches at ATLAS. 

of the LHC to probe the flavour mixing of the heavy neutrinos in left-right symmetric 
models. The dedicated analysis of lepton number violating effects at the LHC requires a 
thorough simulation of the relevant same-sign dilepton background. On the theoretical side, 
it also requires a detailed specification of the lepton number symmetry breaking mechanism. 
Within the minimal left-right symmetric model, lepton number is broken at a high scale, 
generating the heavy Majorana masses of the right-handed neutrinos. The resulting LNV 
effects at the LHC are therefore maximal, with the heavy neutrinos decaying with equal 
probabilities into positive and negative leptons. An overview of the sensitivities of high and 
low energy probes of LFV and LNV processes in this scenario is shown in Figure [161 

Especially with respect to the origin of lepton number violation, the scenario analyzed in 
this work is not unique, and there are theoretical consideration to address some of the issues 
of the minimal left-right symmetric model. Right-handed neutrinos are the messengers 
whose exchange yields neutrino masses through the type-I seesaw mechanism. Similarly, 
heavy scalar triplet exchange induces neutrino masses through the type-II seesaw. It is 
therefore expected that, at some level, the smallness of neutrino masses will make it difficult, 
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if not preclude, to probe the physics of the heavy right-handed neutrinos. First note that 
having TeV-scale right-handed gauge bosons in the minimal SU(2) L <g> SU(2) R <g> \J(1) B _ L 
model discussed in Sec. Ill Al does not, by itself, provide a fully satisfactory picture. For 
example, the gauge couplings in such minimal scheme have no simple SO(10) embedding 
compatible with gauge coupling unification. On the other hand right-handed neutrinos at 
the TeV scale can only be possible through the ad-hoc requirement of tiny Dirac neutrino 
Yukawa couplings in Eq. (EJ). 

These shortcomings can all be naturally evaded by implementing a low scale seesaw mech- 

or linear seesaw [79[ within the SU(2)x 



anism, such as inverse [761 
context 



77j, [78| 



SO 



SU(2) R ® U(1) B _ L 

This would not only justify the lightness of the right-handed neutrinos, 
without the need to invoke unnaturally small Yukawa couplings, but it also achieves a con- 
sistent gauge coupling unification within supersymmetry 87 . 
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